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Abstract 

We study the triton and three-nucleon force at lowest chiral order in pionless effective 
field theory both in the Hamiltonian and Euclidean nuclear lattice formalism. In the case 
of the Euclidean lattice formalism, we derive the exact few-body worldline amplitudes 
corresponding to the standard many-body lattice action. This will be useful for setting 
low-energy coefficients in future nuclear lattice simulations. We work in the Wigner 
SU(4)-symmetric limit where the S'-wave scattering lengths 1 So and 3 Si are equal. By 
comparing with continuum results, we demonstrate for the first time that the nuclear 
lattice formalism can be used to study few-body nucleon systems. 
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1 Introduction 



The ultimate goal of nuclear physics is the derivation of physical observables from the underlying 
theory of the strong interactions, quantum chromodynamics (QCD). At low energies quarks and 
gluons are confined in nucleons and pions and one is forced to work in the non-perturbative 
regime of QCD. A model-independent approach to non-perturbative QCD is provided by lattice 
QCD simulations. However, modern lattice QCD simulations work with lattices not much larger 
than the size of a single nucleon. In the foreseeable future lattice QCD calculations will therefore 
not be suited to obtain direct results in few-body or many-body nuclear physics. 

Alternatively, one can work with the effective degrees of freedom at such energies, i.e., 
the nucleons and pions. Their interactions are summarized in an effective chiral Lagrangian 
which in turn may be utilized, e.g., as the interaction kernel in the Lippmann-Schwinger or 
Faddeev-Yakubovsky equations for few-nucleon systems. The advantage of the effective chiral 
Lagrangian approach with respect to traditional potential models is the clear theoretical con- 
nection to QCD and the possibility to improve the calculations systematically by going to higher 
chiral orders. Over the last few years, effective field theory methods have been applied to few- 
nucleon systems, see e.g. [1,2] and references therein. In particular, in the three-nucleon system 
this framework provided an explanation for universal features such as the Phillips line [3], the 
Efimov effect [4,5], and the Thomas effect [6], which have no obvious explanation in conventional 
potential models [7]. 

Very recently, chiral effective field theory methods on the lattice have been applied to nuclear 
and neutron matter [8-1 1] . The chiral effective Lagrangian is discretized on the lattice and the 
path integral is evaluated non-perturbatively using Monte Carlo sampling. In this work we 
employ pionless effective field theory at lowest chiral order on the lattice to study the triton 
and three-nucleon force. We work in the Wigner SU(4)-symmetric limit where the 1 S'o and the 
3 Si scattering lengths are the same, and the isospin-spin SU(2) xSU(2) symmetry is promoted 
to an SU(4) symmetry. 

This is the first study of the three-nucleon system in the lattice formalism and the motiva- 
tions for this study are twofold. First, it provides a starting point for more realistic studies of 
the triton system as well as other few-nucleon systems using effective field theory and efficient 
Monte Carlo lattice methods. Second, it provides the first measurement of the three-nucleon 
force on the lattice in the Wigner symmetry limit. This provides an important step towards fu- 
ture many-body nuclear simulations with arbitrary numbers of neutrons and protons including 
three-body effects. 

This work is organized as follows. We start in the next section by reviewing continuum 
regularization results for the triton binding energy. This will serve us as a reference for the 
lattice calculations in the subsequent sections. In Sec. 3 we consider the Hamiltonian lattice 
using the Lanczos method, while in Sec. 4 the Euclidean lattice formulation is utilized and the 
path integral is evaluated with Monte Carlo sampling. Our results indicate that many-body 
simulations of dilute nuclear matter in the Wigner symmetry limit should be possible without 
a sign problem, even for unequal numbers of protons and neutrons. This will be discussed 
in Sec. 5, while our conclusions are presented in Sec. 6. A few formulae are deferred to the 
Appendix. 
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2 Continuum regularization 



The nontrivial dependence of the three-body force on the cutoff is a non-perturbative effect. 
There is no finite set of diagrams which produces the ultraviolet divergence. According to the 
Thomas effect [6] a two-body interaction with range R and fixed binding energy B 2 produces in 
the zero range limit a deeply bound three-body state with a binding energy that scales as l/R 2 , 
see [7] for a review. This singular behavior of the three-body system yields a three-body force 
that could be quite a different function of the cutoff scale for different regularization schemes. 
In this section we therefore consider first the usual continuum regularization which will then 
be compared to lattice regularization employed in the following sections. 

In the SU(4)-symmetric limit, where the isospin and spin degrees of freedom can be in- 
terchanged, the lowest order non-relativistic effective Lagrange density including three-body 
interactions reads 



L = ip 




(i) 



where ip is a four-component column vector of the four nucleon states with mass m, 



1> = 



Pi 



(2) 



After summing bubble diagrams in the dinucleon or dimer propagator and renormalizing the 
two-body interaction, we have for the renormalized value of Cq the relation 



^0 — j 



m 



(3) 



where a 2 is the two-body scattering length. Following [12,13], the homogeneous S'-wave bound 
state equation for dimer-nucleon scattering is 



„ \ 2 f t(k, 
t{k,p) = - / dq 

0.2 V 



t{k,q)q 2 



- J 3 4-mE 3 



1 j fp 2 +pq + q 2 — mE 3 
pq \p 2 — pq + q 2 — mE 3 



3mC 2 



(4) 



where A is the cutoff momentum, E 3 is the total energy, and k is determined by 

/4m . I, 
k=^ — (E 3 + B 2 ) , 

with B 2 the two-body binding energy 

B 2 



1 



(5) 



(6) 



In the S'?7(4)-symmetric limit, the 1 S'o and 3 Si two-body sectors are degenerate. Since the 1 S'o 
channel has a nearly zero energy resonance, we take the S'?7(4)-symmetric two-nucleon binding 
energy to be about one-half the deuteron binding, 



Bo = 1 MeV. 



(7) 
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Figure 1: The three-body interaction Dq as a function of the triton energy £3 for various cutoff 
momenta using continuum regularization. The two-nucleon binding energy is set at 1 MeV. 



The coupling Dq is shown in Fig. 1 as a function of £3. In the notation of [12] we have 



rr/ a\ _ hA 2 _ £> A 2 _ £> mA 2 (R . 
1 ' 2mg 2 6mC 2 96vr 2 a 2 ' 1 ' 



with H(A) ~ 1 away from critical values of A, and so -Do scales roughly as A -2 for £"3 much less 
than —1 MeV, as can be seen in Fig. 1. We observe a pole in Dq as a function of £3. For the 
cutoff values considered here the pole occurs rather close to the —1 MeV continuum threshold 
for a dimer plus nucleon. The location of this pole, £f ole , slowly increases in magnitude with the 
cutoff A. For £3 greater than £f° le , there is no way to accommodate the triton as the ground 
state of the three-body system. As we cross E^ ole we can instead switch labels and identify the 
triton with the first excited 5-wave state. If we do this then the value of Dq must switch from 
infinitely repulsive to infinitely attractive, and the new ground state now lies outside the range 
of validity of effective field theory. For very large Aa2 we expect to see many such poles, but 
in our case Aa2 is less than 10 and there is only one pole. 



3 Hamiltonian lattice with Lanczos method 

In this section we consider the triton and three-body force in the Hamiltonian lattice formalism. 
After constructing the lattice Hamiltonian for the three-nucleon system, we use the Lanczos 
method [14] to find the lowest eigenvalues. Since we compute several low eigenvalues, the 
Lanczos method is well-suited to observing the singular behavior of the three-body coupling 
Dq in the triton energy. However, the method requires storing in memory the space of three- 
nucleon states at rest. Thus we cannot probe larger volumes nor any useful sytems with more 
than three particles. 
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Let Oj(n s ) be an annihilation operator for a nucleon with spin-isospin index i at the spatial 
lattice site n s . The lattice Hamiltonian with two-body and three-body interactions correspond- 
ing to the Lagrangian in Eq. (1) is 

3 . 1 r 

H = ~ ^ a l(^s)oj(n s ) - — ^ |aj(n s )a;(n s + Z s ) + aT(ra s )ai(ra s - Z s 

n s ,i n s ,i i s =l,2,3 

+ C ^ a\(n s )a i (n s )a t j (n s )a j (n s ) + D ^ aj(n s )a i (n s )aj.(n s )a j (n s )4(n s )a jfc (n s ) (9) 

n s ,i>j n s ,i>j>k 

with the unit lattice vectors l s in the spatial directions. All parameters are given in lattice 
units, i.e., physical parameters are multiplied by the appropriate power of the lattice spacing 
a. As derived in [11], C is determined by the nucleon-nucleon scattering length a 2 , 

Co = y X » ( 10 ) 

"* 47ra 2 L3 ^ Q_. 

ky^O integer 

where 

n « o 27r/ci o 27r/c2 o 27r/c3 

iL = o — 2 cos 2 cos 2 cos . 11 

k L L L v 1 

One disadvantage of using lattice discretization is that rotation invariance is no longer 

exact. We cannot project out S-wave bound states as in standard continuum regularization. 

Nevertheless we can identify the subspace of three-nucleon states with zero total momentum 

and construct the corresponding Hamiltonian submatrix. We then use the Lanczos method to 

identify the lowest energy eigenvalues and states invariant under the full cubic symmetry group 

of the lattice. 

For lattice spacings ar x = 30,40,50 MeV we consider lattices with length L = 4,5,6,7. 
The dimension of the space is given by L 6 which can be seen as follows: in the rest frame, 
we need only to specify the relative displacements of the three distinguishable particles. This 
is given by two vectors and hence (£ 3 ) 2 = L 6 possibilities. We then use these finite volume 
results to extrapolate to the infinite volume limit. The results are shown in Fig. 2. As with 
the continuum regularization results we see that D scales roughly as A -2 . In contrast with 
continuum regularization results, however, it seems that D is also shifted towards more positive 
values for smaller values of A. There is again a pole singularity not too far below the threshold 
value of —1 MeV and the location of this pole again slowly increases in magnitude with the 
inverse lattice spacing. 



4 Euclidean lattice with Monte Carlo 

In this section we consider the triton and three-body force in the Euclidean lattice formalism. 
We employ the lattice action described in [9] for pure neutron matter and extend it to include 
protons in the Wigner SU(4) limit. Since the formalism was originally designed for many- 
body simulations at fixed chemical potential and arbitrary number of particles, we remove the 
dependence on the chemical potential. We now go through the steps somewhat carefully to 
find the corresponding Euclidean lattice formalism for fixed number of particles. Our rationale 
for this extra care is to preserve the same lattice discretization errors due to the spatial and 
temporal lattice spacings as in [9], so that we can use the few-body results to fix coefficients in 
the many-body simulation. 
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Figure 2: The three-body interaction Dq as a function of the triton energy E$ for various lattice 
spacings using the Hamiltonian lattice formalism. The two-nucleon binding energy is set at 
1 MeV. 



4.1 Free particle 

We start by studying the lattice action for a free single fermion. On the lattice the free 
Hamiltonian for a single fermion species is given by, see Eq. (9), 



H free = — ^2 a\n s )a(n s ) - ^— ^ ^ |^a t (n s )a(ra s + l s ) + a\n s )a(n s - l s ) 

ri s n s Z s =l,2,3 

One can approximate the partition function as a Euclidean path integral 
Zf ree = tr [exp (-pH free )] ~ J DcDc* exp (-S free ) , 



(12) 



(13) 



where we have dropped an irrelevant constant on the left side of Eq. (13). The lattice action 
for the free non-relativistic fermion is given by 



gfree = ^ [ c * + 6) — (l — Qh) C * {n)c{n)} 



h [c*(n)c(n + l s ) + c*(n)c(n — l s ) 



(14) 



with h = at I (2m) and at = a t /a and Grassmann variables c, c*, while the signifies one lattice 
unit in the forward time direction. Note that the relation between the partition function and 
the Euclidean path integral in Eq. (13) is exact for h = 0, so that discretizetion effects in 
temporal direction are minimized. By translating the c field back by one unit in the time 
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direction and using the same notation for simplicity, we can write this in the form 

gfree = ^ [ c *(£) c (fl) _ (1 _ Qh) C * (ft)c(ft ~ 6)] 



-hJ2\ c*(n)c(n + l s -0)+c* (n)c(n - Z s - 0) 



(15) 



We now wish to calculate the path integral by considering explicit fermion worldlines, see 
Fig. 3 for sample worldlines. 
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Figure 3: Shown are sample worldlines on a 1+1 dimensional lattice. The circles denote a hop 
during a time step to a neighboring lattice site. 



Clearly it is not feasible to take into account all possible worldlines even for smaller lattices, 
so we approximate the actual value of the path integral by performing Monte Carlo sampling 
of the worldlines. The contribution of the worldines to the path integral is most conveniently 
calculated in the occupation number basis which we will introduce in the following. 

To this end, let us first simplify the situation by considering the action without spatial 
hopping and fixed lattice site n s , 



S f H ree = c*(n)c(n) - (1 - 6h)c*(n)c(n - 1) 



(16) 



where n is the temporal component of n. So we are considering the physics of a single free 
particle at zero temperature with only one spatial lattice site, i.e., a + 1 dimensional system. 
Next, let us define coherent quantum states 



\ v ) = e" at |0> = (l+V) |0> 
where rj is a Grassmann variable. Then 



a\q) 



ae vat |0) =a{l + r]a ] ) |0) = 77 10) 



(17) 
(18) 



and 



(rj'\ v ) = (0| e^V at |0) = (0| (1 + ffa)(l + r/a f ) |0) 
= 1 + fj'r) = 



(19) 
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Also, in this basis the identity is given by 

J dfjdrie-™ \r]) (r]\ = J dfjdr] (1 - 7777) (1 + a^) |0) (0| (1 + rja) 

= / dfjdrj(— 7777) |0) (0| + / dfjdrj (rjfj) |1) (1| 



= |0)(0| + |1)(1| = 1, (20) 

where |1) is the normalized state with occupation number 1. We now do some reverse engi- 
neering. Let O be the transition operator for the fermion from time slice i to the next time 
slice i + 1. We want the path integral with antiperiodic time boundary 

J DcDc*e W (-St; e ) (21) 

to arise from the trace of an n-fold product O n where n is the number of lattice sites in temporal 
direction 

tr(C^) = J driN^dfjM^-^ ( V n-i\0 -OO-O (r^-i) 

= J driN^dfiN-ie^-™- 1 (m-i\0 J df^^N^e-^-™-* \ m - 2 ) 

x { m _ 2 \ O J dfjN^driN^e-^-^N^ \ m _^ ( m _ 3 \ j^} (^j \ m _^ . ( 22 ) 

Explicit knowledge of the quantum operator O is not necessary, but its matrix elements are 
required to fulfill 

(77,1 O |77i_i) = exp [77,77,-1(1 - 6h)] (23) 

such that the trace of O n , Eq. (22), reproduces the path integral. 

Let us now calculate matrix elements of O in the occupation number basis. For brevity we 
merely present the calculation for (0| O |0), while the results for the remaining matrix elements 
are deferred to the appendix, 

(0| O |0) = (0| J dfj'dr]' e'^' |t/) (t/| O J d^e^ 11 (77) (77I 0) 

dfj'dr]' J df)di)e- m e- r1ri exp [^(l - 6h)] (0| 7/) (r]\ 0) 
dfj'dr]' J dfjdr]e~^ e _f] v exp [fj' 77(1 — 6h)] 
= J dfj'dr)' J dfjdr)(l - 7777) (1 - fj'r)') [1 + fj'r)(l - 6h)} 

= J dfj'dr]' J drjdr](-rjr])(-rj' r]') = 1- (24) 

Next, we consider the case with more than one lattice site and with nearest-neighbor hop- 
ping, 

gfree = ^ [ c *(tf) c (tf) _ (1 _ Qh)c*(n) C (n - 6)] 
n 

-hJ2 [c*(n)c(n + l 8 - 6) + c*(n)c(n - l s - 6)] . (25) 
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Let ?7„_i, rj n be the Grassmann variables at given spatial lattice site (call it A) at time steps 
n — 1 and n. Let 9 n -i, 9 n be the Grassmann variables at a neighboring spatial lattice site (call 
it B) at time steps n — 1 and n. In order to reproduce the path integral we demand 

(r} n , 9 n \ O |?7„_i, = exp [77^77^—1(1 - 6/1)] exp [0 n 0„_i(l - 6/1)] 

x exp [fj n 9 n -ih\ exp [^77^-1/1] . (26) 

We now construct the matrix elements of O in the occupation number basis for sites A and B. 
To this aim, we consider the subspace when there is exactly one particle in either A or B. Let 
|1, 0) be the state when the particle is at A. Let |0, 1) be the state when the particle is at B. 
We find 



<l,0|O|l,0) = J dfj'dr]' J dfjdr] J dB'dB' J d9d9{\ - 7777) (1 - fj'r}') (1 - 99)(1 - 9'9') 
x [1 + f)'r]{\ - Qh)} [1 + 9'9{\ - Qh)] [1 + fj'9h] [l + 9'rjh] rj'fj 
= (1 - Qh). (27) 

Also 

(l,0|O|0,l) = J dfj'dr]' J df)dr) J dB'dB' J d9d9{\ - fjrj) (1 - fj'r}') (1 - 99)(1 - 9'9') 
x [1 + 07(1 - Qh)} [1 + 9'9{l - Qh)] [1 + fj'9h] [l + 9'r}h] rj'9 
= h. (28) 
Clearly by symmetry between sites A and B one has 

(0,l|O|0,l) = (l-6/i) (29) 
(0,l|O|l,0) = h. (30) 

The rules are now clear for a single free fermion worldline: for a hop to a neighboring lattice 
site we get an amplitude h, whereas one obtains (1 — Qh) if the particle stays at the same lattice 
site. 



4.2 Two-body interaction 

We continue by introducing a two-body interaction, the same for every pair of different fermion 
species as required in the Wigner SU{4) limit. We take the two-body contact term of the 
Hamiltonian in Eq. (9) 



H 2 = C ^ a K^s)a i (n s )a t j (n s )a j (n s ) . 



(31) 



n s ,i>j 



Utilizing a Hubbard- Stratonovich transformation with an auxiliary scalar field s(n) eliminates 
the four-fermion term and yields 

Co 



S = J2 [ c *(n)c t {n + 6) - e v - CoatS ^ +£s ^{l - 6/i)c*(r1) C ;(n) 

n,i 

- h [cm4(n + Is) + c*(n)c t (n - L)] + ± £ s\n). 



(32) 



n,l s ,i 
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The strength of the two-body coefficient Co is determined by summing nucleon-nucleon bubble 
diagrams [9]. Again we translate the c field 



S = J2 [c*i{n)ci{n) - e^WW+^l - 6h)c*(n) Ci (n - 6)] 

n,i 

-hY, [<{n) Ci {n + l s - 6) + c*(n) Cl (n - l s - 0)] + ± £ * 2 (*0 



(33) 



and factorize the two-body interaction contribution 

exp [-S] = exp [S free ] 



x 



n 



e 2 



s 2 (n) 



exp 



| e ^ s (n)+fe - l) (1 - 6/i) c *(n)c t (n - 6) 



(34) 



For simplicity we start with the case when there are only two fermion species which we call 
i =|, |. Expanding the exponential of the two-body interaction in Eq. (34) and making use of 



J ds{n) e-¥ 2 ^ ( e v^cw(n) + ^ _ i) = o 



(35) 



we arrive at 



/ 



ds(n) e 2 s2 (") exp 



ds(n) 



i 

e -^ 2 (n) + e -\s\n) ^- CaatS {n)+^ - l)' (1 - 6h) 2 c* (n)c T (n - 6)cJ(n)q(n - 0) 



= v^exp [(e- Coat - 1) (1 - Qhfc\(n)c^(n - 6)cJ(n)q(n - 6)] . 



(36) 



The v^7r is an irrelevant factor which drops out of every physical observable and will thus be 
neglected. The important term is 



exp [(e- Coat - 1) (1 - Qhfc\(n)c^(n - 6)cJ(n)q(n - 6)] 



(37) 



which yields a contribution only if an j particle and a j particle are at the same spatial site at the 
beginning of a time step and neither one hops during the time step. For such a configuration the 
free particle action delivers an amplitude of (1 — Qh) 2 such that the entire amplitude including 
the interaction becomes 



[1 - Qh) 2 + (e~ Caat - 1) (1 - Qh) 2 = e- Coat (l - 6hf 



(38) 



Hence there is an additional factor of e~ Coat due to the two-body interaction. The reason 
for this simple result is that when h — 0, there is no time discretization error in our lattice 
formulation. It exactly matches the Hamiltonian formulation without hopping. 
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We generalize our findings to the case with more than two fermion species. Expansion of 
the two-body interaction term yields 



cxp 
= f + 



^ e V=cW(n)+3p _ l j (1 _ q H) J2 c*(n)c t (n - 6) 

i 

^ e v/=cw(rt)+Sp - l) (l - Qh) c*Mci{n - 6) 

i 

+ ^ e V=?W(n)+3p - l) 2 (1 - Qh) 2 c*(n) Cl (n - 6)c*(n) Cj (n - 6) 

i>j 

+ (V I ^(™)+ £ ^-l) 3 (l-6/i) 3 c*(n) Cj (n-6)c*(n) Cj (n-6)c*(n)c fc (n-6) 

i>j>k 

+ •••, (39) 
where the ellipsis denotes higher orders in the Grassmann variables c*, c. The integral formula 



e _l s 2( f J) + fc v ^Co^ S (fJ) + fc^p 



^e^^v^. (40) 

allows us to perform the integration over the auxiliary field s. For three different fermion 
species, which we denote by % — 1, 2, 3, we obtain after integration over s 

exp [(e- Ci,at - 1) (1 - 6/i) 2 c*(n)ci(n - 0)c* 2 (n)c 2 (n - 6)] 
x exp [(e~ Coat - l) (1 - 6h) 2 c* 2 (n)c 2 (n - 6)c* 3 (n)c 3 (n - 6)] 
x exp [(e' Coat - l) (1 - 6h) 2 c* 3 (n)c 3 (n - 6)c*(n)ci(n - 6)] 

x exp [(e- 3Coat - 3e~ Coat + 2) (1 - 6/i) 3 ct(n)ci(n - 6)c;(n)c 2 (n - 6)c* 3 (n)c 3 (n - 6)] . (41) 

Note that the apparent three-body interaction from the Hubbard-Stratonovich transformation 
is 0(aj), whereas the other terms in Eq. (41) are of order 0(a t ). It is therefore a discretization 
effect which disappears as the temporal lattice spacing goes to zero. 

It is straightforward to see that this term simply produces the correct combinatorial factors 
in the pairwise two-body interactions, and in fact is not a three-body interaction. If there are 
exactly two particles of different species at the same site and neither one hops during the time 
step, we get a contribution analogous to Eq. (38), 

(1 - 6h) 3 + (e- Coat - 1) (1 - 6h) 3 = e- Coat {l - 6h) 3 , (42) 

i.e., we obtain an extra factor of e~ Coat from the two-body interaction. But for three different 
particles at the same site and no hop during the time step the contribution to the amplitude is 

(1 - Qh) 3 + 3 (e- Coat - 1) (1 - Qh) 3 + ( e - 3Coat - 3e- Coat + 2) (1 - Qh) 3 
= e - 3Coat (l-Qh) 3 . (43) 

So we get a factor of e" 3Coat from the three pairwise two-body interactions. Again, the reason 
this is so simple is that when h = 0, there is no time discretization in our lattice formulation. 
Moreover, from these considerations it should be clear that for an arbitrary number of fermion 
species we get a factor e~ c ° at for each pairwise interaction. 
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4.3 Three-body interaction 

In the present investigation we choose the lattice three-body interaction in such a way that for 
exactly three particles of different kind at the same site and no hop during the time step the 
contribution to the amplitude reads e~ D ° at (l — 6h) 3 . We assume that the lattice action with 
Hubbard-Stratonovich fields has been designed appropriately to coincide with this definition. 
Its explicit form is not relevant here. 

4.4 Results 

After having discussed the possible contributions to the worldlines of three distinct fermions 
on the lattice, we present now the results of the simulation. Let |0, 0, 0) be the state with three 
nucleons, each of different kind and each with zero momentum. We employ path integral Monte 
Carlo with worldlines to compute lattice approximations for 

G{(3) = (0, 0, 0| exp [-/3#] |0,0,0) , (44) 

making use of the fact that the state |0, 0, 0) has a non-vanishing overlap with the three-body 
ground state of the Hamiltonian H . By measuring the exponential tail of G{(3) for large (3, we 
have a measurement for the triton ground state energy E 3 . For pairs of spatial and temporal 
lattice spacings, 

a' 1 = 30 MeV, a" 1 = 9 MeV; (45) 
a' 1 = 40 MeV, a" 1 = 16 MeV; (46) 
a' 1 = 50 MeV, a" 1 = 25 MeV; (47) 
a' 1 = 60 MeV, a" 1 = 36 MeV; (48) 

we compute Do as a function of the triton energy E 3 . The results are shown in Fig. 4. 

We have stayed away from the sensitive region E 3 ps — 1 MeV. Since there are many dimer 
plus nucleon continuum states near the energy threshold —1 MeV, it is difficult to extract E 3 
accurately. It requires lattices with larger temporal extensions and significantly more CPU time 
to study this critical region. For E 3 < —2 MeV, however, the simulation is not difficult and 
the finite volume effects are small. For a" 1 = 30 MeV and 40 MeV we used a lattice volume 
of L 3 = 6 3 ; for cr 1 = 50 MeV we used L 3 = 7 3 ; and for cr 1 = 60 MeV we used L 3 = 8 3 . 
The error bars shown are stochastic error estimates, determined by the data from four different 
processors with completely independent runs. The finite volume errors are significantly smaller 
than the stochastic error estimates shown. We note the clear similarity between the Euclidean 
lattice results and Hamiltonian lattice results. 

5 Lattice action positivity 

Pure neutron matter has been studied on the lattice using effective field theory both with 
pions [8] and without pions [9-11]. For cold dilute neutron matter at momentum scales below 
the pion mass, pionless effective field theory should provide an adequate description of the low- 
energy physics. One nice feature of the lowest-order pionless theory for pure neutron matter 
is that it can be implemented on the lattice with a positive semi-definite action by means of a 
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Figure 4: The three-body interaction Dq as a function of the triton energy E$ for various 
lattice spacings using the Euclidean lattice formalism. The two-nucleon binding energy is set 
at 1 MeV. 



Hubbard-Stratonovich transformation. This feature is very important for performing efficient 
Monte Carlo simulations. 

For cold dilute nuclear matter with a small proton fraction, one expects that pionless ef- 
fective field theory also describes the relevant low-energy physics. However, in this case a 
three-nucleon force is required for consistent renormalization. Until recently it has been an 
open question whether or not this three-nucleon force would spoil positivity of the lattice ac- 
tion. It was shown in [15] that in the Wigner SU(4)-symmetric limit, the lattice action would 
remain positive so long as the three-nucleon force was not too strong and the four-nucleon force 
was not too repulsive. However, it was not known if the parameters of the real world satisfy 
these conditions. In this section we show that these conditions are indeed satisfied for typical 
lattice spacings relevant for cold dilute nuclear matter. 

Consider a Euclidean lattice action with two-, three- and four-body interactions of the form 

h(^(n)) 2 + c^(n)f + c 4 (^(n)) 4 ] . (49) 

Ft 

In [15] it was shown that a Wigner SU(4)-symmetric nuclear lattice simulation is possible with 
unequal numbers of protons and neutrons without a sign problem if and only if the matrix 

"10 2c 2 

2c 2 6c 3 
_2c 2 6c 3 12c 2 2 + 24c 4 

is positive semi-definite, with the added condition that if c 2 = then c 4 = 0. The determinant 
of this matrix is lGdj — 36C3 + 48c 2 c 4 . The relations between c 2 and C3 (in lattice units) and 



(50) 
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a' 1 (MeV) 


a; 1 (MeV) 


h 


C (MeV" 2 ) 


D (MeV" 5 ) 


c 2 


c 3 


16c^ - 36c 2 ; 


30 


9 


0.0532 


-2.66 x 10~ 4 


4.52 x 10~ 9 


0.28 


-0.016 


0.35 


40 


16 


0.0532 


-1.83 x lO" 4 


2.22 x 10~ 9 


0.25 


-0.023 


0.23 


50 


25 


0.0532 


-1.39 x lO" 4 


1.09 x 10~ 9 


0.23 


-0.026 


0.18 


60 


36 


0.0532 


-1.13 x 10~ 4 


5.64 x 10~ 10 


0.22 


-0.027 


0.15 



Table 1: Shown are the values for c 2 , c 3 and the determinant 16d] — 36c§ for given values of the 
lattice spacings a and a t . 



our parameters Co and D (in physical units) are 

c 2 = i(l-6/i) 2 (e- c ^-l), (51) 
c 3 = 1(1 - Qh) 3 (e- Doa ~ 6at - l) . (52) 

While the three-nucleon force is required for consistent renormalization at lowest order, the 
four-nucleon force is an irrelevant operator and can be neglected for sufficiently small lattice 
spacings [16]. If we assume the four-nucleon force to be zero, then the determinant is 16c|— 36c§. 

For an SU(4)-symmetric triton binding energy of 8 MeV, we find the results shown in 
Table 1. We see that for lattice spacings in this range the determinant 16^ — 36c 2 remains 
positive. This implies that many-body simulations of cold dilute nuclear matter in the Wigner 
symmetry limit should be possible without a sign problem, even for unequal numbers of protons 
and neutrons. 

As noted above, we have ignored the irrelevant operator producing an SU(4)-symmetric 
four-nucleon force. If we do include a small four-nucleon force, this positivity result should 
most likely remain intact. However, this will be addressed explicitly in a future study of the 
four-nucleon system [17]. 

6 Conclusions 

In this work we have introduced two novel approaches to few-body nuclear lattice simulations. 
We worked in the Wigner SU(4)-symmetric limit where the S'-wave scattering lengths are equal. 
First, in the Hamiltonian lattice formalism with the Lanczos method the three-body interaction 
strength was calculated as a function of the triton binding energy for various lattice spacings, 
while keeping the two-nucleon binding energy fixed at 1 MeV. Our findings indicate that the 
coupling scales roughly — as expected — as the lattice spacing squared. We furthermore observe 
a pole singularity in the coupling slightly above the threshold value of 1 MeV in the triton 
binding energy. This is in good agreement with continuum regularization results. 

The second framework applied is the Euclidean lattice formalism with Monte Carlo sampling 
of the path integral. We observe a clear similarity with the Hamiltonian lattice results for 
triton binding energies greater than 2 MeV. For binding energies closer to the threshold energy 
of 1 MeV it is difficult to extract the binding energy accurately with this method due to the 
presence of many dimer plus nucleon continuum states in this energy region. 

Although we have restricted ourselves to three nucleons and two- and three-nucleon forces 
in the present investigation, our results suggest that the Euclidean lattice method can be 
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generalized to a larger number of nucleons and more complicated forces amongst them. It 
sets the stage for further studies in the few-body sector, such as the inclusion of a four-body 
force which can be easily included in this formalism or effects due to breaking of Wigner 
symmetry [17]. The importance of higher chiral orders in the action can be studied as well. 

Moreover, this work provides a first measurement of the three-nucleon force on the lattice 
in the Wigner symmetry limit and is an important step towards future many-body simulations 
with arbitrary numbers of nucleons including three-body effects. 
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In this appendix, we present the remaining matrix elements of O in the occupation number 
basis not given in the main text. For a lattice with one spatial lattice site and one fermion we 
obtain 



A Matrix elements of O 
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